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A VERSATILE MOSSBAUER SPECTROMETER AND ITS APPLICATIONS 
IN VIBRATION MEASUREMENT 

By Jag J. Singh and Lona M. Howser 
Langley Research Center 

SUMMARY 

The present study was undertaken in order to design, construct, and test a bench 
model of a Mossbauer spectrometer efficiently operable in the scattering mode. Such a 
spectrometer has been developed and used to perform experiments intended to demonstrate 
the applicability of Mossbauer-effect spectroscopy to important problems in structural 
mechanics without any special specimen preparation. One such experiment which dealt 
with the measurement of low frequency (<10 Hz) and low amplitude (<0.254 mm) periodic 
vibrations of austenitic steel specimens is described in detail. From these investigations, 
it is concluded that using the present spectrometer, measurements can be made on a 
variety of structural materials under environmental conditions without any special speci- 
men preparation. The spectrometer system has been constructed almost entirely from 
commercially available components. 


INTRODUCTION 

In order to extend the applicability of the Mossbauer effect to nondestructive mea- 
surements in surface stress studies, structural vibration analysis, and liquid-solid inter- 
actions, it is essential to use reflection (backscattering) technique instead of the standard 
transmission method. Backscattering technique, however, has not yet been fully developed 
due to the special detector configuration requirements for a large solid angle at the surface 
under examination. 

Several attempts at backscattering technique have been recently reported (refs. 1 
to 4). These studies are based on the detection of either the reemitted 14.4 keV gamma 
rays or the 6.3 keV conversion X-rays or the conversion electrons. Although the Fe^- 
backscattering technique should be most efficient when the low energy conversion elec- 
trons are detected, it is suitable for very thin surfaces (<3000 A) only. Even for thin- 
surface studies, this technique requires that the specimen be kept in high vacuum and is 
further hampered by the limited penetrating power of the low-energy conversion electrons 
through the detector windows. A new Fe57 -radiation detector, which can subtend an angle 
as large as -it steradians at the sample surface in both the transmission and the reflec- 


tion geometries, has been designed and developed. The detector is relatively rugged in 
construction since it involves the detection of the reemitted 14.4 keV gamma rays and the 
6.3 keV conversion X-rays. 

The detector has been used in the measurement of low frequency (<10 Hz) and low 
amplitude (<0.254 mm) cyclic vibrations of austenitic steel specimens. When a Co57_pt 
source is moving with constant-acceleration motion, the Mossbauer spectrum from a sta- 
tionary scatterer is characterized by a single peak centered at a certain source velocity 
p Q . If now the scatterer is subjected to low frequency, low amplitude periodic motion, 
which has no phase relationship with the source motion, conditions of Mossbauer resonance 
between the emitter and the absorber atoms can occur at any source velocity in the range 
of p Q ± V max where Vmax represents the maximum scatterer velocity. This situation 
leads to a characteristically modified Mossbauer resonance peak from a moving scatterer. 
A comparison between the spectra from the stationary and the moving scatterers can then 
be used to infer the particulars of the scatterer motion. (Even though, by its very nature, 
the Mossbauer effect gives information about the relative motion of the source -absorber 
atoms, there have not been many applications of Mossbauer spectrometry to vibration anal- 
ysis as such. The few reported techniques are of the contacting type, applicable over lim- 
ited velocity ranges and depend upon the use of standard reference absorbers. (See refs. 5 
and 6.) One such study (ref. 6) typically involved obtaining velocity calibration absorption 
spectra of various source-absorber combinations. The source-absorber combinations 
were subdivided into three cases: resonance, partial resonance, and off-resonance. All 
of these combinations involved attaching the Mossbauer source to the system under test and 
using standard-reference stationary absorbers in transmission geometry. Such measure- 
ments, besides involving source mounting adaptation, require electronic gating techniques 
to limit counting intervals to selected velocity ranges. In order to eliminate these require- 
ments, a new, noncontacting, technique of adapting Mossbauer -effect spectroscopy to vibra- 
tion measurements in materials containing Mossbauer nuclei has been developed. (See 
ref. 7 for the list of these nuclei.) This technique does not require any special specimen 
preparation, is applicable over a wide velocity range (submillimeter to several 
centimeters/second) and lends itself easily to determining the type of specimen-motion 
waveform. It has been applied to analyze the vibration of stainless steel specimens.) 
Results of several experiments involving different types of specimen motion and computer 
programs developed for analyzing the complex spectra, observed when the source and the 
scatterer are in random simultaneous motion, will be described. 

SYMBOLS 

A amplitude of Lorentzian peak 


2 



I 


I(x) observed counts in the channel centered at source velocity x (constant 

acceleration) 

i index 

* = V max /AV 

p Q position of the Lorentzian peak on the source-velocity axis corresponding to 

zero scatterer velocity 

S. - AV 
1 dV/dt 


t time 

V(t) scatter velocity at time t 

V m ax maximum value of the scatterer velocity 

x source velocity 

r full width of Lorentzian peak at half maximum intensity 

co angular velocity 

Subscripts: 

o,i indices 


RADIATION DETECTOR AND ITS PERFORMANCE 

The radiation detector consists of a 1.524- by 5.08- by 5.08-cm Nal(Tl) crystal with 
a central hole 1.524 cm in diameter. A sodium iodide scintillation crystal was chosen as 
the detector because of its high efficiency for low energy gamma (or X-) rays and easy 
availability in various geometrical shapes and sizes. The crystal encapsulation incorpo- 
rates front and back beryllium windows 0.0127 cm thick and 3.81 cm in diameter, that per- 
mit easy penetration by the conversion X-rays and the reemitted 14.4 keV gamma rays 
into the sensitive region of the Nal(Tl) crystal. At the two ends the crystal is coupled, 
via quartz light guides, to two matched photomultipliers. Figure 1 shows a schematic 
drawing of the detector assembly. The electronic gains of the two photomultipliers are 
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equalized and the signals at the respective collectors are added and amplified before 
passing through a single channel analyzer (SCA). The limits of the SCA are set to accept 
all pulses in the range of 5 — 15 keV. Figure 2 shows a simplified diagram of the elec- 
tronic circuit. The collimator assembly shown in figure 2(b) keeps the incident photon 
beam divergence at the scatterer to an angle less than 2.5°. This low divergence reduces 
the source velocity dispersion at the scatterer to less than 0.1 percent. The radiation 
counter, in combination with a sensitive electromechanical transducer and drive system 
producing motion of the source with constant acceleration through a range of velocities, 
has been used to study the Mossbauer spectra in both the transmission and reflection 
geometries. 1 Figures 3 and 4 show the transmission and the reflection spectra, respec- 
tively, obtained with a 0.8 mg/cm2, Fe^-enriched iron target. Figure 5 shows a typical 
backscattered Mossbauer spectrum obtained with a 0.0762 -cm thick, nonmagnetic stainless 
steel (type 347) scatterer. The Mossbauer source in all these studies, was a 10-millicurie 
Co^ source in platinum matrix. The absorber or scatterer in each case was positioned 
to subtend an angle of it steradian at the detector (see fig. 2(b)). This large solid angle, 
besides reducing the experimental run time, improves the signal-to-noise ratio in the 
spectrum since part of the background is independent of the solid angle subtended by the 
detector at the scatterer, that is, scattering at the collimator edges and surrounding 
shield. From the spectra shown in figures 3 and 4, it is clear that the present detector 
configuration is capable of producing equally good transmission and reflection Mossbauer 
spectra. The observed resolution in the present study is largely limited by the Nal(Tl) 
detector resolution for the 5 to 15 keV electromagnetic radiation. 

For an efficient use of the Mossbauer spectrometer for various types of studies 
indicated earlier, a computer program for analyzing experimental spectra has been devel- 
oped. In this program, the experimental spectra are least-squares fitted to a curve made 
up of a parabola and a series of Lorentzian peaks superimposed on it. (See ref. 8.) For 
those cases where the Lorentzian peaks were smeared — such as those obtained when the 
scatterer had a sinusoidal or constant acceleration motion — a modified version of the 
program, described later, was used. 


APPLICATIONS 

The Mossbauer spectrometer has been used in a number of studies in structural 
mechanics. One such study is described below. This study involved measurement of 
the bulk motion of steel specimens subjected to low frequency (<10 Hz) and low amplitude 
(0.254 mm) periodically varying velocities. 2 Three types of periodic motions were 

1-The velocity scale of the spectrometer was calibrated using a standard soft iron 
transmission Mossbauer spectrum. 

^These limits were determined mainly by the source characteristics (strength and 
recoilless fraction) and the scatterer-detector solid-angle considerations in the present 
study. 
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imposed on the specimens using a standard laboratory waveform generator.. The speci- 
mens were 0.0762- by 2.54- by 10.16-cm pieces of austenitic stainless steel (type 347), 
mounted vertically at the end of a horizontal shaft of an electromagnetic vibrator ener- 
gized by the waveform generator. It was assumed that the specimen motion followed the 
generator waveform. Figure 6 shows a schematic drawing of the experimental setup used 
in this study. The scatterer was subjected to the following types of motion: 

Sinusoidal motion: V(t) = V max sin wt 
Constant-velocity motion: V(t) = ±V m ax 
Constant-acceleration motion: V(t) = ±V max T ft 
where f = Constant acceleration 

and +v max > v (t) > -Vmax 

During all these experiments the source was moving with a constant-acceleration motion 
in the velocity range of ±6(1.274) mm/sec. The various combinations of source -scatterer 
motions are illustrated in figure 7. 

Scatterer Moving With Sinusoidal Motion 

A backscattered Mossbauer spectrum from a stationary nonmagnetic stainless steel 
(type 347) specimen was first taken and is shown in figure 5. However, when the specimen 
is subjected to a sinusoidal motion, which has random phase relationship with the source 
motion, the resonance peak position p Q can fall anywhere between p Q ± V max where 
V max corresponds to the maximum value of the scatterer velocity along the direction of 
Co57_gource motion. Thus, by measuring the spread in the resonance peak position, the 
maximum velocity of the scatterer along the source -scatterer direction can be computed. 
An independent measurement of the frequency of the scatterer motion should then permit a 
determination of its displacement (amplitude). Figure 8 shows a typical Mossbauer spec- 
trum obtained when both the source and the scatterer were in motion with random phase 
relationship. It is clear from this figure that the scatterer spends maximum time at a 
velocity close to the maximum value, as expected for a sinusoidal motion. (Fig. 8 also 
suggests the possibility of using this technique to determine the exact waveform of the 
scatterer motion.) Figure 8 and other similar spectra, shown in figure 9, corresponding 
to different scatterer velocities, were computer analyzed by least-squares fitting the 
observed data to an expression of the form 

2n 

I(x) = (ax 2 + bx + c) + ^> Sj (1) 

IS (m - x) 2 + (r/2) 2 
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for various 


where p. = + i XsifS S; represents relative dwell times / - ^ V \ 

10 n 1 \dV/dty 

constant -velocity intervals (AV) in the range (-V max — +V max ), and n = max . Fig- 

' ' A V 


ure 10 illustrates the graphical technique of calculating Sj. The computer program, 
used in this analysis, was a modified version of the basic program described in refer- 
ence 8. The modifications, along with the solution of a sample case, are described in the 
appendix. A comparison between the values of ±V max , computed from the Mossbauer 
spectra and from the measured values of the frequency and amplitude of the scatterer 
motion is shown in figure 11 and summarized in table I. (The velocity was calculated 
assuming the scatterer motion to be sinusoidal in response to the input sine wave.) The 
agreement is quite good. 


Scatterer Moving With Constant- Velocity Motion 

The nonmagnetic stainless steel (type 347) scatterer was subjected to a constant - 
velocityS motion in the range ±3(1.274) mm/sec and backscattered Mossbauer spectra 
were taken, in the same manner as described previously, for a constant-acceleration 
source motion in the range ±6(1.274) mm/sec. Figure 12 shows typical spectra for this 
type of source and scatterer motions. Notice the appearance of two additional peaks on 
either side of the zero-scatterer-motion peak. Spectra shown in figure 12 and others cor- 
responding to different values of scatterer velocities were least- squares fitted 4 to an 
expression of the following form: 


I(x) = (ax2 + bx + 






it! (Pi - xy + (r/2r 


( 2 ) 


where 


p l = p o,l “ p o " v max 


and 


p 3 = p o,2 ~ p o + V max 

3 The constant-velocity calibration was obtained using Mossbauer spectroscopic 
technique. (Direct measurement.) 

4 The computer program used in this analysis was essentially the same as that 
described in reference 8, except for the fact that there were three independent resonance 
peaks — instead of the usual one (nonmagnetic) or six (magnetic) for uniphase steel 
spectra. The outside peaks were least-squares fitted using a half-width value corre- 
sponding to a stationary scatterer whereas the central peak required a slightly larger 
value due to its velocity dispersion. 
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The scatterer velocity amplitude is obtained from the position of the two extreme peaks 
p Q i and p Q g as follows: 

Po,l = p o " V max 
p o,2 = p o + V max 

where p Q is the peak position for a stationary scatterer. (It is also the position of the 
central peak.) 

The central peaks in figure 12 arise from the fact that the scatterer velocity goes 
through zero in the (±V max — ^Vmax) transition in a finite (i.e., nonzero) time. A com- 
parison between the Mossbauer-computed and directly measured scatterer velocities is 
shown in figure 13 and is also summarized in table I. The agreement between the two sets 
of values is quite good. 

Scatterer Moving With Constant-Acceleration Motion 

If the nonmagnetic steel scatterer is subjected to a constant-acceleration motion 
simultaneously with the constant-acceleration source motion, and if there is no phase cor- 
relation between the two motions, the resultant backscattered Mossbauer spectrum is 
expected to be of the following form: 

2n 

I(x) = (ax2 + bx + c) + Y — (3) 

i=l (Pi - x ) 2 + (r/2 ) 2 


where 




Equation (3) is identical to equation (1) if Sj were the same for all velocity intervals 
AV. The Mossbauer spectra for this combination of source and scatterer motions can be 
computer analyzed using the same program as that for sinusoidal motion described in the 
appendix. 

Figure 14 shows an experimentally observed spectrum with the scatterer moving 
with constant-acceleration motion. A theoretically predicted spectrum under these condi- 
tions is shown by a solid line in this figure. It is obvious that the experiment and the 
theory do not agree in this case. The experimental spectrum indicates that the scatterer 
spends more time at or near extreme velocities than elsewhere. An examination of the 
waveform-generator output showed that it was distorted (see insert in fig. 15). The theo- 
retical spectrum was then computed using an expression of the type given in equation (1) 
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with the Sj values approximate to the scatterer-motion waveform in the present case. 
The S i values used are shown in figure 15. Figure 16 shows a comparison between the 
experimental and the theoretical spectra under these conditions. The agreement is rea- 
sonably good. 

This exercise has demonstrated that it is possible to extract information about the 
scatterer-motion waveform from the Mossbauer spectra observed when both the source 
and the scatterer are in random relative motion. From the Sj values required to fit 
the observed spectra, one can reconstruct the scatterer-motion waveform. Conversely, 
Mossbauer-effect spectroscopy can serve as the basis of control circuitry for maintaining 
a desired scatterer waveform over long periods of time. 

Figure 17 shows a theoretical modification superimposed on a nonresonant 
Mossbauer spectrum in order to indicate how ±V m ax °f the scatterer can be calculated 
from constant-acceleration scatterer-motion spectra. In this figure, the source and the 
scatterer motions had the following values: 

Source velocity (constant acceleration) = ±6(1.274) mm/sec 

Scatterer velocity (constant acceleration) = ±2(1.274) mm/sec 

Notice the raised section in the central part of the parabola. The peak positions corre- 
sponding to P 0 (beginning) and p o(end) give a measure of v max b Y the following 
relation: 

T 7 - _ ^o(end) _ P o (beginning) 

•max - 2 

CONCLUDING REMARKS 

A radiation-counter configuration suitable for reflection Mossbauer spectroscopy 
has been described. An application of the spectrometer, involving measurement of peri- 
odic vibrations of steel specimens has been discussed. The measured and the computed 
values are in agreement within the expected limits of errors. The errors on the computed 
values range from 0.7 percent to 3.0 percent in the case of sinusoidal motion and 0.4 per- 
cent to 1.6 percent in the case of constant-velocity motion. The major source of error in 
the computed values lies in the statistical accuracy of determining the extreme scatterer - 
velocity points in the Mossbauer spectra. This error can be reduced by using a stronger 
source (i.e., 100 milli curies Co^ or stronger) which will also have a further desirable 
effect of reducing the experimental time correspondingly. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., June 21, 1972. 
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APPENDIX 


MODIFICATION OF THE COMPUTER PROGRAM DESCRIBED IN 
REFERENCE 8 TO ALLOW FOR THE EFFECTS OF 
SCATTERER MOTION AND A SAMPLE CASE 

Computer Program 

. The basic computer program D3290 developed at Langley Research Center to ana- 
lyze a complex Mossbauer spectrum from a stationary scatterer has been described in 
detail in reference 8. This program has been slightly modified to include the effects of 
scatterer motion in the present study. The modified program D3291 least-squares fits 
the observed spectrum to an expression of the form given in the text in equation (1), which 
describes a parabola plus a series of gradually displaced Lorentzian peaks. The instruc- 
tions for the use of program D3291 are the same as for the basic program D3290 except 
for a few additions. These additions and modifications are described below. 

The program D3291 is run as if there was a single peak. The input NP = 1 and 
the dimensions have been reduced to allow for only one peak. Three new inputs are 
required arid two inputs are redefined. An additional variable is printed with the output 
and two are redefined. These are listed below: 




New Input 

FORTRAN variable 

Dimension 

Description 

CGAMO 

(i) 

Peak width at half maximum intensity, set equal to 
the half- width value for a stationary scatterer. 
(D3290 gives the value for a stationary scatterer.) 

NDT 


Number of subdivisions of the velocity cycle 

PO 

(i) 

Initial estimate for p Q - V max 

P20 

(i) 

Initial estimate for p 0 + V max 

RT 

(1300) 

Relative dwell times, Sj 



New Output 

Heading 


Description 

P 


Position of (p Q - V max ) on the velocity axis 

P2 


Position of (p Q + V max ) on the velocity axis 

CGAM 


Peak width at half maximum intensity 
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APPENDIX - Continued 


The results given in this report were obtained by dividing the scatterer motion into 
1248 subdivisions. This number was selected to ensure that the corresponding V max 
had saturated to its true value within the limits of accuracy of the present system. The 
program can be altered to find an optimum number of subdivisions for a given experi- 
mental system. 

Program D3291 requires approximately 60 000 octal locations of core storage. A 
typical case using 1248 displaced Lorentzian peaks runs in 400 seconds. 

The listing for program D3291 is as follows: 

PROGRAM 03291' ( INPUT=201 , OUTPUT = 201 , TAPE 5= INPUT ,T A PE6=HUTPUT , 

1 SC FILE *TAPE7=SCFILE) 

COMMON A { 1 1 , A0(1 ) , AREA! 1 ) , ARRAY { 1025, 7) ,P( 7 ,8 ) ,C(7, 1 ) ,CGAM( 1) , 

2 CCAM0( i ) , ERFOP (7J » P( 1) » P0( 1) ,P2 ( U t P2Q( 1 ) , P ( 1025) t RT (1300 J , 

2 PS(650)»X(1025)»Y{ 1025) 

DIMENSION C0NPY< 1025 ) ,PAR (1025) ,YYPAR ( 2050 ) 

DIMENSION PFM{ 3>,XM(3),YM(2) 

GOUT VALENCE ( PPM ( 1 ) , S AG ) , ( PRM ( 2 ) , EC ) , ( PR V ( 3 ) , CD ) 

EQUIVALENCE ( C C M p Y ( 1) » A R R A Y ( 1*1 ) ) , ( P AR (1) , 4 FR AY ( 1 , 2 > I 
•EQUIVALENCE ( YYPAP ( 1 ) * ARRAY ( 1 , 3 ) ) 

DIMENSION IPIVCT ( 2 1 ) » I NDEX( 21*2) 

EQUIVALENCE ( R ( 1 ) , I P I V0T ( 1 ) ) , ( R ( 22 ) , INDEX! 1,1)) 

DIMENSION V( 7,7) 

EQUIVALENCE (V (1,1), ARRAY (1,4)) 

NAMELIST /NAM1 / AO , BO , CG AMC , EP RC R , I ERR , IF L AG, I PR IN T , NOT , NP , PO , P2C , 

1 PT,SAn,vELR 

DATA XM , Y m /27H SOURCE VELOCITY (MM/SEC), 19H COUNTS PER CHANNEL/ 

NPC = 7 

CALL CALCOMP 
CALL LEROY 
OCIN = l 

2 READ (5, KAMI) 

TF (EOF, 5) 3,4 
3 CALL OALPLT (0,0,999) 

STnp 

4 CONTINUE 

CC-AM( 1 )=CGAMC'( 1) 

C INTERPOLATE, GET MORE DWELL TIMES 
N1 = 40 

5 CONTINUE 

DQ 6 T = 1 , N 1 

6 RS<n = RT(I) 

DO 7 I = 2 , N 1 
K = 2 *1 -l 

RT(K-l) =( PS(I-l) * RS(I))/2. 

7 PT(K)=RS( I ) 

A'1 = K 

IF ( N 1 .LT.NDT) GQ TO 5 

NDT =K 

ANDT=K 

I F ( I FLAG.NE. 1 ) GO TO 27 
READ (5,9) I D , N C 
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APPENDIX - Continued 


9 FHPMAT (5XtI5fT10) 

IF (NO). 10,1?, 12 

10 *'C = -NP 

CPM =-1 

11 FORMAT ( 8(F4*0,F6.0) ) 

12 PEAD (5,11) ( X( I ) ,Y( I ) ,1=1 ,NO) 

C SCAM AP P AY AND INTER POL AT E FOR RAC POINTS 
J=1 

DP 18 T = 2 , NO 

IF ( Y ( I ).EG*0. ) GH TO 16 
IF (J.EC.l) GO TO 18 
FAP= ( Y ( I ) — Y ( K ) ) / J 
..t = J-l 

DO 10 N as 1 , J 

13 Y(K+N)=Y<K-1+N)+ FAC 
J = 1 

. or TO 1 8 

16 IF ( J * F C * 1 ) K=I-1 
J = J+1 

1 8 'CONTINUE 
PEW IND 7 

WRITE (7) ( Y{ I ) * I ~1 , NO ) 

rn=Y( l ) 

ANO = N 0- 1 

xi=xm 

SUM = VELP *2* 

DP 21 T = 1 , N 0 

X U ) = (VELP-SIJV* (X ( I )-Xl )/ ANO) *CCN 
21 CONTINUE 

27 *TEP=0 

WPTTF (6,2700) I D , I ER R 
2700 FPPmat ( * 1 C A $ E NO • * 1 5 , 5X* I EP R = * I 3 ) 

T F ( T EP P • EQ* 1 ) GO TO 32 
pEWTvd 7 

PFAD { 7 ) ( V( I ) ,1 = 1 ,N0 ) 

r n TP (32,28,29), IEPR 

28 PH 2310 1=1, NO 

Y(I ) = Y { I )+ SORT(Y(I ) ) 

2810 CONTINUE 
GO TP 32 

29 00 2910 I = 1 , A 0 

YU ) = Y ( I )- SQR T ( Y ( I ) > 

2910 CONTINUE 

3? DP 40 T = 1 , NO 
XSO=X( I )**2 
A P p A Y ( T , 1 ) = X S 0 
A PP AY ( I » 2 ) = X( I) 

A PP AY ( T , 3 ) = 1 • 

FX= SAC*XSQ + R0*X ( I ) + CO 
PIT )- ( Y ( I ) - FX) 

40 CONTINUE 
80 NP3=3*NP +3 
N 1 = N P 3 + 1 

C INITIALIZE D LCT ROUTINE AND SET ORIGIN 
m,m = 4 

qh no P = 1 1 N P 
SAVE=CGA YC (M)**2 /4*0 
SAVET = AO( N 1 ) * CGAN"0(N)/2.0 
A2= 2.0 *AC(m) 

P 6 = (P?P( V)-PC(M) )/ ANDT 



APPENDIX - Continued 


DO 100 T — 1 ? NO 

PM X = pom - xm 

sa=o.o 

sp=o*o 

SG= 0 • 0 

SP2=0.0 

DP 90 J=1,NDT 

PNU^= PM> +J*P6 

DENG^PMJM **2 + SAVE 

DSQ= QENCN*DENO 

S A= S A + P T ( J ) /DENG V 

SG= SO - PT(J)*SAVE1 /DSQ 

AJ=J/ANCT 

PAPT=PT( J)*A2*FNUM/DSQ 
SP=SP -PAFT*< l.-AJ) 

SD2 = .S 0 2- PART*AJ 
90 r0NT T NUE 

ARRAY ( T ,VM) = SA 
ARP AY ( T » Y M + l ) =SP 
ARRAY ( I , m|v !+2) =S P 2 
FX= A D ( M ) *APRAY(I,MM» 

R ( T >-P ( I )-FX 
100 CONTINUE 

M M = V M + 3 

110 continue 

C ARRAY TRANSPOSE * ARRAY 
DO 120 K= 1 f NP 3 
HP 138 ^ - 1 » N P 3 
P(Kf M )=0*0 
on 115 T = 1 tMG 

PCKfM)=ARPAY(I,K)*ARRAYCItM] + 6(K,M) 

115 rnNTINUE 
118 CONTINUE 
120 CONTINUF 
C ARRAY TRANSPOSE * R 
on 1 2 q ^ - 1 » N P 3 

r (k ) = o.o 

on 125 ^ = 1 » NO 

125 C{K ) = AFQAY<M,K )*R m + C(K) 

129 CONTINUE 

r ALL MATINV (B f NP3,C, l f DfcT, IPIVCT, I MCE X * 7 

ITEP=ITER+L 

DO- 14? 1=1,3 

TES T 1 = ABS<C( I >/PRM( I) ) 

IF (TEST1.GT*EPP0R( I ) ) GO TO 300 
142 CONTINUF 
N M= 4 

DO 145 w=l,NP 

TEST1 = ABS(C(MY)/AO(M) ) 

JF (TEST1.GT .ERROR C^M)) GO TO 800 
TEST 1 = ABS(C(MN + 1)/PG(Y)) 

TF (TEST1.GT*E C R0R( W M+1) ) GO TO €00 
TEST 1 = rBS(C(H^+? >/ P2J(M)) 

IF { TES T 1 *GT •€ RROR C MM+2 I ) GO TO £00 
145 CONTINUE 
165 ww =4 

SAP=SA°+ C(l) - 
B0= BP +C (2 ) 

rn=cn + cm 


f I SCALE ) 
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APPENDIX — Continued 


00 168 V=1,NP 
A ( M ) = Cl v IP)+ A C ( N ) 

P(M)= f (NN+1) + P0( N) 

P? ( M ) =C ( VN +2 ) ♦ P2CHM) 

W|V|= MM + 3 

168 rnNTTNUF 

on 170 1=1, NO 
qn=Y(i] 

170 rnviPY( T ) = 0.0 
on 174 V=1,NP 

P6= (P?r { M)- PC(N))/ ANDT 
S AVE = 0 GAM N J **2/4.0 
On 172 1*1 i NO 

PMX = P ( M >— Xtll 
on 171 J= 1 » NO T 

0 = P6*J 

op A K =FT(J)*A(N)/( ( PMX+D >**2+SAVE) 

OQMPY ( T ) = CCMPY (I ) + PEAK 

171 TONTTNUF 

172 CONTINUE 
174 CONTINUE 

on 176 T = 1 » N P 

p ar 1 1 ) = san*x( n**2 + ao*x( i ) +cc 

fPMPY(I)= COMPYl I ) + PAR(I) 

176 p { T )=R ( T)-C,rNPY(II 
STO=0.0 
on 178 T = 1 1 NO 
178 STD= D ( T )**2+ STD 
ANN*NO-N P3 
STD* S TO / ANN 
DP 181 T=1 »NP3 

on 130 J = 1 * NP 3 
Vd ,J)=P< I » J )*STD 

180 GPNT! NUE 

181 CONTINUE 

DO 183 T=1»NP3 
183 \M T tl ) = SOPT( V l I» I) ) 

STO = SOP T ( STD ) 

P 1 = VP 1.8** 3/3 . 

P 3= VE L B**2/2. 

AR 1 = S AO* P 1 + BO*P 3 + CP*VELB 
AR2= -SA0*P1 + B0*P3 - C 0*VEL3 
ARE Al= A R 1 - A R 2 
AR E A? =0 . 

WRITE (6,185) SAP, V( 1,1) , RO » V ( 2 , 1 ) , CP , V< 3 , 1 ) » I TEP 
185 FORMAT (4 8HOO0EFFIC IENTS OF PARABOLA Y= SA*X**2 + B*X *■ C/ 

1 SA=*E15.8, *( T*E9.3»*)*,5X*B = *E15.8, * ( T*E9 . 3 , * )* , 5X ,*C=*E 15. 8 » 

3 1F + ,19X, 1H— , 33X t 1H— , 33 X » 1H— /*0N0 . OF I TERATI ONS=* 16 ) 

DP 200 N = 1 »NP 

EVALUATE INTEGRALS TO GET AREAS FOR PHI2 

SAVE = C0 AN ( M> **2/4.0 
»6= (P2(M)-P(M)) /ANDT 

OP 191 J = 1 »ND T 
01= A(«)*PT(J) 

A J = J 

02 = SAVE +P(M)**2+ 2 -0* P ( M ) * AJ*P6* (AJ*P6)**2 
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APPENDIX — Continued 


ac = c 2 

BC=— 2 . 0 * ( A J# P6+ P(M) ) 

CC = 1.0 

0 = 4.0 *AC*CC — BC **2 
IF CO ) 186,940,167 

186 FAC = 1 2 . *CC + PC 
FAC1= SOFT (- C) 

FAC 2= (FAC - FAC 1 ) /{FAC + FAC1) 

FAC =A|_rC ( F AC 2 ) 

API = FAC / FAC1 

FAC. = -12. *CC + BC 

F AC 2= (FAC — F AC1 ) / (FAC + FAC 1 ) 

FAC = ALOG (FAC2) 

AR2 = FAC/ F AC 1 
AREA(M)={ AR1- AP 2 ) *C1 
ACE A 2= AFEA(M) + AFEA2 
GO TO 189 

187 FAF1= SORT (0) 

FAC =( 12.0*CC + PC ) 

FAC 2= A T Afv2(FAC,FACl ) 

AR1 = 2.C * FAC. 2 / FACl 
FAC =(-12.0*CC +BC ) 

FAC2=ATAN2(FAC,FAC1 ) 

AP 2 = 2.0* FAC 2/ FACl 
AP F A ( M ) = { AR 1— AR2) *C1 
ARE A2= AP F A( M ) + A R E A 2 

189 CONTINUE 
191 CONTINUE 

PHI 2= APS ( AREA2 /AREA1 ) 

ATS=P(M) 

WRITE (6,190) M 

190 FORMAT ( # OPAP A RE TER S FOR PEAK*I3) 

WRITE (6,195) A I S 

195 FORMAT (*0IS=*E15 .8 ) 

M $ a 3 * M 

V| n S+4 ) *C • 

WRITE (6,205) A(M) , V( MS+1 ) ,P( M ) , V(MS+2) ,P2( N) , V(M$+3) ,C GAM ( M ) , 

1 V ( M S + 4 ) 

WRTTF (6,204) 

200 WRITE (6,203) AP£A(M) 

203 FORMAT ( *OARfc A=*E 13.6) 

204 FORMAT ( 1H + , 1 6 X, 1H- , 3 IX , 1H- , 32X , 1H- , 34X , 1H- ) 

205 FORMAT ( *0A=* E 1 3 . 6 , * ( T* E9 . 3 , * ) * , 5X* P = *E 1 3. 6 , *(T’f t E9.3,*)*5X’i tD ?=* 

1 E13.6 , *( T*E9 .3,*)*, 5X*CGAM=*E 1 3 . 6 , *{ T*E9. 3, * )* ) 

C EVALUATE THESE PARAMETERS IF ONLY CAE PEAK 
YMAX=SAC*P(M)**2 * RO*P(m) + CC 
PHI 1 = A ( N ) /YMAX 
POP =P H I 1/PHI 2 

WRITE (6,206) PH 1 1 , PH I 2 , POP 

206 FORMAT (*O p HI 1=*E15.8»5X*PHI2=*E15.8»5X*PHI 1/PHT2**615.8) 

210 IF (IPRINT. EQ.O) GC TO 380 

WRITF (6,230) 

230 FCIPMAT (*OCHANNEL N0.*3X*X MM/S E C* 14X*Y* I4X *C n "Pi*TEO Y*10X 
1 *RFSIDLALS*10X*PAP ABOLA*) 

DO 315 1=1, NO 

WRITE (6,312) I , X ( I ) , Y( I ) ,CGMPY ( I ) ,P (I ) ,PAR(T) 

312 FORMAT ( I 6 ,5E2 0. 8 ) 

315 CONTINUE 

WPTTE (6,370) STD 
370 FORMAT (*OSTD=*E 15.6) 
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APPENDIX — Continued 


380 K = 1 

IF ( TEP's.l^E*l ) GC TO 2 
C COMPUTE M I N T MJMS AND MAXIMJMS 
XP0=l 2 • 

TAIL ASrALt (X tXPGfNO»Kf 10* ) 

YPG = 1. 0 # 
f- UM=2*NP 
on 400 T= 1 ♦ NO 
400 Y Y P A P (I ) = Y ( I ) 

N*i =MP + i 

on 410 1= Ml t MUM 
• 410 YYPAP{ T )=PAR( I-NO) 

• TALL ASCALE ( YYPAR » YPG#NUM , K , 10 • ) 

NP 1=N0+ 1 
NP2=NP4-2 

Y(MP1 ) = Y YP AP ( NUM + 1 ) 

Y(MP? )=YYPAR ( NUM+2 ) 

PAO (NPl ) = Y ( M P 1 ) 

PAP(NP2)=Y<NP2) 
rpM py ( Nj P 1 ) =Y ( N P 1 ) 
rrMDY(MP2)=Y(MP2) 

C HP AW X AXIS 
XQV=10. 

X T I C = 1 • 

TAIL AXES (O.,0o0. ,XPG f X(NPll f X < NP 2 > . XT I C , XO V , XM t # 1 5 ,-27 ) 

C DRAW Y AXIS 
Y D V = 1 0. 

YT T f = 1 • 

CALL AXES (0. , 0 . , 90 . , YP G , COMP Y ( N PI ) ,CCMPY (NP2 ) * YTI C v Y DV, YM , .1 5 y 19 ) 
C 

C PLOT CURVE 
r 

CALL PL PT ( X t Y » NO ) 

IF (ITFP.GE.30 ) GO TO 600 

TAIL ITNPLT (XfCOMPY f MOtK f O,OtOfO) 

CALL LTNPLT ( X , PAP f NO , K , 0 , 0 t 0 , 0 ) 

C ESTABLISH A MEW REFERENCE POINT FOP THE NEXT GRAPH 
600 CALI CALPLT (l4«*0.'-3) 

GP TP 2 
800 M M=4 

IF (ITEP.GE.3i)) GO TO 1003 

CP 900 N = 1 1 N P 

AP ( M) = r ( YM) + AC ( M) 

PP( M) =C ( YY + 1 ) + P 0 ( R ) 

B2 n (M) =c (NM+2) + P2CIM) 

995 wm=mv+3 

900 CONTINUE 
920 SA^=SAP+ C ( 1) 

F Ps BP+ C ( 2 ) 

CO =rc+ C ( 3) 

PO TP 32 

940 WRITE ( 6 » 950 ) A * D * CG AM 

950 FORMAT (*OWHEN EVALUATING THE INTEGRALS C ECUALS 0., THIS CANNOT 
1 BE t A»PfCGAM ARRAYS FOLLOW * / ( 7 E 1 6 . 6 ) ) 

GO TP 2 

1000 OP 1020 ! = 1 » N 0 
1020 PAR (I )=Y( I) 

CP TP • 380 
FND 
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APPENDIX - Continued 


Sample Case 

The input and output for a sample case are printed below. In this case, the source 
was moving with a constant acceleration in a velocity range of ±6(1.274) mm/sec. At the 
same time, but with no phase correlation with the source motion, the scatterer was under- 
going a sinusoidal motion in the range of ±V max = (1.16 ± 0.01) mm/sec. The source 
and the scatterer motions were collinear. The computer printout does not include a 
velocity calibration factor of 1.274, 

The input for the sample case is as follows: 


*NAM1 IFLAG=l, I PR INT=0 » IERR=1 » $AO=50., B0=-5., EPR0P=2*l.E-5, l.E-2, 

3*1. E-5, 

NP=1 , A 0= . 3 * P0=-1. 1 » P 20= .6 * CGAM0=.2929* VFLB=6., NDT=1248, 

P T=. 95 , . 35, . 285, .25,. 225, .205,. 19, .175, .17, .165, .16, .155, .151, .15*. 149,. 146, 
.145, .141, .14, .139, 

.139,. 14,. 141,. 145,. 146, .149,. 15,. 151,. 15 5,. 16,. 165,. 17,. 175,. 19,. 205,. 225, 
.25, .285 ,.35, .95, 


$ 


2 

4 

37200 

3 

-1004 

37214 

4 

3721 1 

5 

37557 

6 

36875 

7 

37155 

8 

37092 

9 

36589 

10 

36978 

11 

37096 

12 

37016 

13 

36710 

14 

36787 

15 

36930 

16 

37214 

1 7 

36369 

18 

36770 

19 

36778 

20 

37274 

21 

36937 

22 

36988 

23 

37161 

24 

36931 

25 

37211 

26 

3701 7 

27 

37057 

28 

36934 

29 

36865 

30 

36541 

31 

36891 

32 

36836 

33 

36626 

34 

36389 

35 

36500 

36 

36500 

37 

36554 

38 

36653 

39 

36996 

40 

36709 

41 

36724 

42 

36 419 

43 

36944 

44 

37000 

45 

36842 

46 

36624 

47 

37000 

48 

36694 

49 

36700 

50 

36779 

51 

36500 

52 

36493 

53 

36712 

54 

37000 

55 

36923 

56 

36712 

57 

36731 

58 

36513 

59 

36587 

60 

36053 

61 

36415 

62 

36834 

63 

36557 

64 

3.65 71 

65 

36494 

66 

36663 

67 

36659 

68 

36676 

69 

36777 

70 

36596 

71 

36979 

72 

36711 

73 

36600 

74 

3662 0 

75 

36710 

76 

36637 

77 

36600 

78 

36630 

79 

36662 

80 

36472 

81 

36887 

82 

36641 

83 

36339 

84 

36812 

85 

36601 

86 

35996 

87 

36115 

88 

36843 

89 

36858 

90 

36 48 5 

91 

36585 

92 

0 

93 

36255 

94 

36546 

95 

36535 

96 

36449 

97 

36320 

98 

36796 

99 

3663C 

100 

36679 

101 

36435 

102 

36080 

103 

36467 

104 

36551 

105 

36223 

106 

36498 

107 

36377 

108 

36356 

109 

36655 

110 

36302 

1 11 

36316 

112 

36272 

113 

36245 

114 

36126 

115 

36258 

116 

36199 

1 17 

36546 

118 

36470 

119 

36414 

120 

36489 

121 

36498 

122 

36129 

123 

36225 

124 

36293 

125 

36170 

126 

36368 

127 

36466 

128 

35973 

129 

35989 

130 

36344 

131 

36275 

122 

36192 

133 

36297 

134 

36463 

135 

36478 

136 

36069 

137 

.36171 

138 

36540 

139 

36445 

140 

36007 

141 

36419 

142 

35928 

143 

36283 

144 

36160 

145 

36208 

146 

3 5 94 8 

147 

36012 

148 

36473 

149 

36469 

150 

36259 

151 

36082 

152 

36212 

153 

36028 

154 

36 34 7 

155 

36092 

156 

35915 

157 

36583 

158 

35375 

159 

36388 

160 

36030 

161 

35929 

162 

36203 

163 

36114 

164 

36125 

165 

36050 

166 

35952 

167 

35918 

168 

36088 

169 

36317 

170 

36269 

171 

35642 

172 

35897 

173 

36471 

174 

35749 

175 

35825 

176 

36121 

177 

36398 

178 

36085 

179 

36301 

180 

35858 

181 

36254 

182 

35665 

183 

36249 

184 

36291 

185 

35651 

186 

35853 

187 

36000 

188 

35991 

189 

35986 

190 

36032 

191 

35854 

192 

35684 

193 

36025 

194 

36045 

195 

36166 

196 

36356 

197 

36031 

198 

35885 

199 

35800 

200 

35895 

201 

35781 

202 

3502 8 

203 

36094 

204 

35968 

205 

36108 

206 

36244 

207 

36043 

208 

36109 

209 

35880 

210 

35 74 1 

2! 1 

35705 

212 

36052 

213 

35832 

214 

3571 8 

215 

35850 

216 

35916 

217 

36000 

218 

3600 C 

219 

35719 

220 

35844 

221 

35634 

222 

35542 

223 

35585 

224 

35704 

225 

36055 

226 

35923 

227 

35796 

22 8 

35985 

229 

35707 

230 

35631 

231 

35822 

232 

35603 

233 

35697 

234 

35917 

23 5 

3561 5 • 

236 

3 5 90 1 

237 

35759 

238 

35933 

239 

35773 

240 

36014 

241 

35496 

242 

35976 

243 

35775 

244 

35794 

245 

35825 

246 

35447 

247 

35662 

248 

36054 

249 

35488 

250 

35895 

251 

35972 

252 

35670 

253 

35510 

254 

35850 

255 

36148 

256 

35793 

257 

35789 

258 

35 56 3 

269 

35672 

26 0 

35833 

261 

35569 

262 

35784 

263 

35952 

264 

35890 

265 

35837 

266 

35366 

267 

0 

268 

35565 

269 

36012 

270 

35655 

271 

35584 

272 

35954 

273 

0 

274 

35704 

275 

35441 

276 

35666 

277 

35654 

278 

35824 

279 

0 

280 

35572 

281 

35686 

282 

3583 1 

283 

35859 

284 

3561 8 

285 

35577 

286 

35816 

287 

3591 1 

288 

35757 

289 

35740 

290 

35674 

291 

35436 

292 

35685 

293 

35600 

294 

35727 

295 

36046 

296 

35479 

297 

36030 

298 

35773 

2 G 9 

35416 

300 

35404 

301 

35529 

30 2 

3551 1 

303 

35494 

304 

35759 

305 

35852 

306 

3552 3 

307 

35829 

3C8 

35762 

309 

35671 

310 

35541 

31 1 

35789 

312 

35239 

313 

35886 
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314 

35598 

31 5 

35524 

316 

APPENDIX 
35713 317 35756 

— Continued 
318 35426 319 

35305 

320 

3556*2 

321 

35584 

322 

35438 

323 

3 5524 

324 

35153 

325 

35663 

326 

35670 

327 

35587 

328 

35336 

329 

35565 

330 

35795 

331 

35581 

332 

35546 

333 

35452 

334 

35733 

335 

35677 

336 

35571 

337 

35723 

33 8 

35702 

339 

35668 

340 

35631 

341 

35543 

342 

35463 

343 

35420 

344 

35424 

345 

35292 

346 

35392 

347 

35499 

348 

35689 

349 

35597 

350 

35386 

351 

35483 

352 

35681 

353 

35442 

354 

3547 1 

355 

35663 

356 

35471 

357 

35672 

358 

35715 

359 

35238 

360 

35279 

361 

35561 

362 

35453 

363 

35492 

364 

35710 

365 

35763 

366 

35527 

367 

35678 

368 

35539 

369 

35598 

370 

35528 

371 

353 )R 

372 

35466 

373 

35225 

374 

35574 

375 

35641 

376 

35091 

377 

35320 

378 

35616 

379 

35080 

380 

35398 

381 

35643 

382 

35761 

383 

35770 

384 

35774 

385 

35897 

386 

35458 

38 7 

35479 

388 

35913 

339 

35706 

390 

35756 

391 

35882 

392 

35473 

393 

35 577 

394 

35955 

39 5 

35961 

396 

35588 

397 

35919 

398 

36012 

399 

35658 

400 

359 79 

401 

35924 

402 

36067 

40 3 

36037 

404 

35884 

405 

36117 

406 

35925 

407 

35880 

408 

36093 

409 

36102 

410 

36294 

41 1 

36001 

412 

36392 

413 

36273 

414 

35799 

415 

36037 

416 

36233 

417 

35993 

41 3 

35993 

419 

^6 063 

420 

36247 

421 

36104 

422 

36484 

423 

36264 

424 

36227 

425 

36070 

426 

36 107 

427 

36390 

423 

36Q9 8 

4 29 

35581 

430 

36206 

431 

36030 

432 

36186 

433 

35886 

434 

36 09 9 

435 

36245 

436 

36119 

4 37 

35548 

438 

36210 

439 

36268 

440 

36247 

441 

36158 

442 

36 114 

44 3 

35985 

444 

36277 

4 45 

36218 

446 

36199 

447 

35950 

448 

36108 

449 

36119 

453 

36030 

451 

36217 

452 

35792 

453 

35898 

454 

36311 

455 

35997 

456 

36119 

457 

359L4 

458 

36095 

459 

3 6 1 9 C 

460 

36069 

461 

36179 

462 

35917 

463 

35995 

464 

35926 

465 

35718 

466 

3620 2 

46 7 

3 58 66 

468 

36021 

4 69 

36119 

470 

35987 

471 

36057 

472 

36457 

473 

36050 

474 

35812 

475 

36154 

476 

36027 

477 

36197 

478 

35871 

479 

35738 

480 

35846 

481 

35901 

482 

36062 

48 3 

36287 

484 

36238 

485 

35961 

486 

35883 

487 

35846 

488 

35719 

489 

36067 

49Q 

36277 

491 

35561 

492 

35 90 C 

493 

36203 

494 

35772 

495 

35777 

496 

35925 

497 

35898 

493 

35813 

4Q9 

36165 

500 

35919 

501 

36033 

502 

35885 

503 

36279 

504 

35993 

505 

35891 

506 

36015 

507 

361 C3 

508 

35668 

509 

35784 

510 

35548 

511 

36113 

512 

35911 

513 

35598 

514 

36245 

515 

35887 

516 

36158 

517 

35814 

518 

36062 

519 

35994 

520 

35966 

521 

35847 

522 

36362 

523 

36 C 85 

524 

35870 

525 

36043 

526 

36298 

527 

35991 

528 

36084 

529 

35884 

530 

35898 

531 

36109 

532 

36100 

533 

35857 

534 

36114 

535 

36060 

536 

36296 

537 

36369 

538 

36044 

539 

36282 

540 

36042 

541 

36C86 

542 

36125 

543 

36168 

544 

36163 

545 

36059 

546 

36375 

*547 

36172 

548 

36118 

549 

35586 

550 

35829 

551 

36217 

552 

36011 

553 

36084 

554 

35979 

555 

3593 1 

556 

36035 

557 

35788 

558 

35808 

559 

35600 

560 

35837 

561 

36066 

56 2 

3583 1 

563 

36116 

564 

35710 

5 65 

35570 

566 

35640 

567 

35922 

568 

35609 

569 

35667 

570 

35944 

571 

35597 

572 

354? C 

5 73 

35663 

574 

35317 

575 

35804 

576 

35560 

577 

35600 

578 

35436 

579 

35287 

580 

35305 

581 

35542 

582 

35608 

583 

35576 

584 

35609 

585 

35328 

586 

35715 

537 

35334 

588 

35 79 8 

539 

35249 

590 

35182 

591 

35486 

592 

35444 

593 

35648 

594 

35382 

59 5 

35434 

596 

35153 

597 

34981 

59 8 

35276 

599 

35631 

600 

35646 

601 

35355 

602 

35248 

60 3 

35455 

604 

35338 

605 

35615 

606 

35468 

607 

35563 

608 

35374 

609 

35486 

610 

35687 

61 1 

35192 

612 

35752 

613 

35258 

614 

34962 

615 

35573 

616 

35704 

617 

35126 

618 

35329 

619 

35975 

620 

35386 

621 

35506 

622 

35178 

623 

35216 

624 

35612 

625 

35192 

626 

35359 

627 

35322 

628 

35429 

6 29 

35572 

630 

35460 

631 

35428 

632 

35393 

633 

35253 

634 

35566 

635 

353 14 

636 

35530 

637 

35359 

638 

35597 

639 

35494 

640 

35551 

641 

35587 

642 

35175 

643 

35465 

644 

35388 

645 

35594 

646 

35562 

647 

35560 

648 

35739 

649 

35341 

650 

3502 8 

651 

^5647 

652 

35651 

6 53 

35526 

654 

35710 

655 

35374 

656 

35535 

657 

35303 

653 

35380 

659 

35765 

660 

3545 1 

6 61 

35633 

662 

35570 

663 

35695 

664 

35542 

665 

35322 

666 

3566 1 

667 

3541 8 

668 

35437 

669 

35645 

670 

35562 

671 

35962 

672 

35715 

673 

35478 

674 

35526 

675 

35603 

676 

35639 

677 

35436 

678 

35873 

679 

35822 

680 

35546 

68 1 

35479 

632 

35 799 

683 

36071 

684 

35567 

685 

35692 

686 

35585 

687 

35613 

688 

35524 

689 

35653 

690 

35607 

691 

35390 

692 

35370 

653 

35875 

694 

35629 

695 

35417 

696 

35576 

697 

35819 

698 

35689 

69 9 

35758 

700 

35610 

701 

35989 

702 

35435 

703 

35447 

704 

35710 

705 

35446 

706 

35523 

70 7 

35530 

7C8 

35659 

7 09 

35347 

710 

35438 

711 

35536 

712 

35584 

713 

35786 

714 

35541 

715 

35686 

716 

35731 

717 

35673 

718 

35569 

719 

35902 

720 

35785 

721 

35817 

722 

35478 

723 

35376 

724 

35540 

7 25 

35548 

726 

35754 

727 

35515 

728 

35386 

729 

35501 

730 

36091 

731 

35574 

732 

35695 

733 

35448 

734 

35454 

735 

35816 

736 

35565 

737 

35965 

738 

3557 1 

739 

35434 

740 

35579 

741 

36056 

742 

35526 

743 

35852 

744 

35483 

745 

36074 

746 

35997 

747 

35456 

748 

35518 

7 49 

35576 

750 

35455 

751 

35677 

752 

35504 

753 

35665 

754 

35800 

755 

35919 

756 

35920 

757 

35472 

758 

35833 

759 

35939 

760 

35528 

761 

35740 

762 

35794 

763 

35573 

764 

35898 

765 

35782 

766 

35758 

767 

35956 

768 

36404 

769 

35840 

770 

35798 

771 

35655 

772 

35669 

773 

35864 

774 

35971 

775 

35653 

776 

35692 

77 7 

35902 

778 

35291 

779 

35728 

780 

35516 

781 

35698 

782 

35672 

783 

35912 

784 

35941 

785 

35621 

786 

35960 

78 7 

35914 

788 

35516 

789 

35936 

790 

35904 

791 

36293 

792 

35705 

793 

35723 

794 

35447 

795 

35906 

796 

36073 

757 

35840 

798 

35832 

799 

35885 

800 

35689 

801 

35614 

802 

35823 

803 

35758 

804 

35745 

305 

36020 

806 

36039 

807 

36007 

808 

36209 

809 

35764 
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APPENDIX — Continued 


810 35812 811 35936 812 36233 813 35828 

818 36092 819 36055 820 36289 821 35651 

826 36094 827 36576 828 35896 829 36099 

834 35902 835 36062 836 36210 837 36186 

842 36004 843 36120 844 36340 845 36226 

850 36264 851 36025 852 36384 853 36167 

858 36125 859 36180 860 35935 861 36000 

866 35893 867 36578 868 36542 869 35955 

874 36253 875 36207 876 36703 877 35984 

882 36198 883 36257 884 35990 885 36410 

890 36155 891 36121 892 36180 893 36014 

898 36601 899 36435 900 36187 901 36332 

906 36433 907 36282 908 36335 909 360C7 

914 36557 915 36732 916 36505 917 36191 

922 36526 923 36411 524 36603 925 36274 

930 36617 931 36427 932 36673 933 36509 

938 36607 939 36534 940 36476 941 36511 

946 36387 947 36782 948 36692 949 36463 

954 36687 955 36492 956 36483 957 36951 

962 36793 963 36768 964 36586 965 36664 

970 36894 971 36612 572 36839 973 36558 

978 36899 979 36772 980 36806 981 36722 

986 36852 987 36536 988 37099 989 36599 

994 36891 995 37164 996 37081 997 37283 

1002 366641003 369911004 369711005 36812 


814 35798 815 35723 816 35897 817 35966 
822 35976 823 36106 824 36359 825 35807 
830 36000 831 36275 832 35977 833 36064 
838 36144 839 36299 840 35943 841 35964 
846 35857 847 36166 848 36265 849 36221 
854 36208 855 36107 856 36174 857 35952 
862 35916 863 36140 864 36018 865 36402 
870 36064 871 36324 872 36061 873 36082 
878 36167 879 35912 880 36159 881 36441 
886 35938 887 36358 888 35997 889 36123 
894 36139 895 36243 896 36582 897 36466 
902 36339 903 36359 904 36518 905 36487 
910 36217 911 36503 912 36412 913 36366 
918 36571 919 36695 920 36333 921 36243 
926 36398 927 36955 928 36542 929 36500 
934 36545 935 36236 936 36582 937 36866 
942 36498 943 36464 944 36603 945 36556 
950 36621 951 36435 952 36432 953 36717 
958 36991 959 36490 960 36599 961 36873 
966 36936 967 36508 968 36599 969 37250 
974 37056 975 36678 976 36874 977 36707 
982 36903 983 36682 984 36767 985 36884 
990 36899 991 36742 992 36637 993 36778 
998 37335 999 367931000 368221001 37020 
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The output for the sample case is as follows: 

CASE NC. 4 I F p P = 1 

COEFFICIENTS OF PAPARnLA Y= SA*X**2 + B*X * C 

SA= 4.64546504E+01 (T6.579E-01 ) B = -5 . I 7309 394E +00 (• I . 699E +00 ) C= 3. 532?0824E + 04m. 210E+01 ) 

NO. OF ITERATIONS* e 

PARAMETERS FOR °E AK 1 
IS=-1 • 1583308 1E+0Q 

A- 3.036201E-OUT8.663E-C3) P=- i . 1 5 833 IE +0 C t T1 . 646E-02 ) P2= 6.092466E-01 (T 1. 447E-02 > CGAK= 2 .929000E-0 1 (TO. ) 

ARE A= 6 • 09036 3E+00 

PHI1= 8.57916 2 ?3E— 06 PHI2= 3 . 793432 85E-03 PHI1/PHI2= 2 .26 1 58 273E -03 

The results are shown in figure 8. 
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TABLE I.- COMPARISON BETWEEN THE CALCULATED AND MOSSBAUER 
COMPUTED VALUES OF THE SCATTERER VELOCITIES 


Sinusoidal scatterer velocity 

Constant scatterer velocity 

Directly calculated values, 
mm/sec 
(a) 

Computed values, 
mm/sec 

Directly measured values, 
mm/sec 
(b) 

Computed values, 
mm/sec 

±(0.467 ± 0.010) 

±(0.487 ± 0.015) 

±(0.500 ± 0.005) 

±(0.497 ± 0.008) 

±(0.701 ± 0.010) 

±(0.719 ± 0.012) 

±(1.000 ± 0.010) 

±(1.003 ± 0.011) 

±(0.939 ± 0.010) 

±(0.957 ± 0.012) 

±(1.500 ± 0.015) 

±(1.501 ± 0.015) 

±(1.158 ± 0.010) 

±(1.126 ± 0.015) 

±(2.000 ± 0.020) 

±(1.979 ± 0.014) 

±(1.401 ± 0.010) 

±(1.382 ± 0.014) 

±(2.500 ± 0.025) 

±(2.501 ± 0.011) 

±(1.635 ± 0.010) 

±(1.630 ± 0.016) 

±(3.000 ± 0.030) 

±(3.023 ± 0.020) 

±(1.864 ± 0.010) 

±(1.820 ± 0.017) 



±(2.098 ± 0.010) 

±(2.076 ± 0.021) 



±(2.339 ± 0.010) 

±(2.377 ± 0.019) 




a The errors on directly calculated values include the combined effect of errors on a 0 and u> in the 
expression, V max = a Q w. (a 0 and co represent the amplitude and the angular velocity of the scatterer motion, 
respectively.) 


k The errors on these values represent the calibration accuracy of the constant-velocity values of the scat- 
terer drive system (namely, 1 percent of the velocity scale setting). 



PHOTOMULTIPLIER 







(b) Schematic diagram of the signal conditioning circuit for 
Mossbauer spectrometer. 

Figure 2.- Simplified diagram of the electronic circuit. 
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" 6 ' 5 ' 4 " 3 '-2 -I 0 I 2 3 4 5 6 X1.274 

source velocity, mm/sec 

Figure 3.- A typical transmission spectrum from a thin (0.8 mg/cm2) Fe 57 -enriched iron 
absorber with the detector described in this paper. (Notice the characteristic scatter- 
ing peaks instead of the transmission dips.) 
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Figure 4.- A typical backscattering spectrum from a thin (0.8 mg/cm2) Fe^ -enriched iron 
scatterer with the detector described in this paper. (Notice the pronounced scattering 
peaks characteristic of scattering geometry.) 
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Time — ► 

Waveform Scatterer Motion Waveform 

Sinusoidal scatterer motion. 



>urce Motion Waveform Scatterer Motion Waveform 

(b) Constant-velocity scatterer motion. 
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(c) Constant-acceleration scatterer motion. 
Figure 7.- Summary of source-scatterer motions. 
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Figure 8.- Typical Mossbauer spectrum obtained when the source and the scatterer were in random 
relative motion. The source was moving with a constant acceleration when the stainless steel 
(type 347) scatterer was subjected to a sinusoidal motion. 






Figure 10.- Numerical calculation of Sj for sinusoidal scatter er motion. 
(The number of velocity subdivisions shown in this figure is arbitrary 
and is for the purpose of illustration only.) 
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COUNTS /CHANNEL 



Figure 12.- Typical backscattered Mossbauer spectra from a 0.762 -cm thick 
stainless steel (type 347) scatterer moving with constant velocity. The 
Co57 

-source had a constant-acceleration motion of ±6(1.274) mm/sec. 
(There was no phase correlation between the source and scatterer 
motions.) 




DIRECTLY MEASURED SCATTERER VELOCITY, mm/sec 


Figure 13.- Comparison between the computed and directly measured values of 
scatterer velocity. The scatterer motion was of constant-velocity type and 
had random phase with respect to constant acceleration motion of the 
source. 
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80100 
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Figure 14.- Experimental Mossbauer spectrum observed from a stainless steel (type 347) scatterer moving 
with constant-acceleration motion. The solid line shows a theoretically predicted spectrum under these 
conditions. 



(Arbitrary Units) 




Scatterer Velocity 

Figure 15.- Relative dwell times Si for the experimental waveform of 
the scatterer motipn (see insert) used to obtain the experimental spec- 
trum shown in figure 14. 
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COUNT PER CHANNEL 




GO 

00 


SOURCE VELOCITY LIMITS ■ + 6(1.274) mm/sec 
SCATTERER VELOCITY LIMITS = + 2 ( 1.274) mm/sec 

THE SOLID LINE REPRESENTS THE EFFECT OF THE SCATTERER MOTION SUPERIMPOSED ON BASE-LINE PARABOLA. 



Figure 17.- Computed backscattered Mossbauer spectrum from a nonmagnetic steel scatterer subjected to a 
true constant-acceleration motion at the same time when the source was also moving with constant- 
acceleration motion. 
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